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We perform a likelihood analysis of the recent results on the anisotropy of Cosmic Microwave 
Background Radiation from the BOOMERanG and DASI experiments to show that they single out 
an effective number of neutrinos in good agreement with standard Big Bang Nucleosynthesis. We 
also consider degenerate Big Bang Nucleosynthesis to provide new bounds on effective relativistic 
degrees of freedom N u and, in particular, on neutrino chemical potential £ a . When including 
Supernova la data we find, at 2<r, N v < 7 and -0.01 < £ e < 0.22, |£ M , T | < 2.6. 
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I. INTRODUCTION 

New results on Cosmic Microwave Background Ra- 
diation (CMBR) anisotropy from BOOMERanG |],|, 
MAXIMA H[|, and DASI @] experiments represent an 
extraordinary confirmation of our present understanding 
of some of the key features of the evolution of our uni- 
verse. The clean evidence for the first acoustic peak of 
temperature anisotropies for CMBR around I ~ 200 Jl]-||] 
strongly supports the scenario of a post-inflationary flat 
universe. On the other hand new results on the second 
and third peak confirm the adiabatic inflationary model 
prediction of acoustic oscillations in the primeval plasma 
driven by gravity, and shed new light on how energy den- 
sity is distributed among several components. This is a 
crucial piece of information which affects many indepen- 
dent cosmological observables, so it is reasonable to ex- 
pect that it will be possible in the next years to have a 
rather clear picture of which cosmological model is actu- 
ally realized in our universe. 

In this respect BOOMERanG and MAXIMA first data 
release [[Ip] already stimulated a wide number of stud- 
ies p|]7|-|15||, aimed to constrain the values of the energy 
density parameters normalized to the critical density, Ot, , 
VL m and Q\, due to baryons, dark matter and an effective 
cosmological constant, respectively. In particular many 
authors have addressed the issue of a tension between 
the determination of Q, b h 2 from CMBR data and Stan- 
dard Big Bang Nucleosynthesis (SBBN) |Jl|,|T|,|l§|o) . 
In fact, the finding of a suppressed second peak in the 
CMBR anisotropy resulted in a rather large value for 
this parameter, fl b h 2 = 0.032±g;gg| at 68% CL 0, while 
the experimental data on primordial 4 iJe and D abun- 
dances, prefer smaller values, £l b h 2 = 0.0I9+q'qq2j -R< e f- 
@, and fl b h 2 = 0.020±0.002, Ref. @ (see also 
at 95% CL. These estimates are obtained assuming three 
standard neutrino degrees of freedom. 

New experimental data from BOOMERanG has re- 



fined the data at larger multipoles, and now single 
out a smaller value for the baryonic fraction, fl b h 2 — 
0.021+q qq3 This is mainly due to an increase in the 
analyzed dataset (roughly by a factor 8) and a better 
understanding of the experimental beam, calibration and 
pointing. The new analysis leads to a slightly increased 
amplitude for the second peak (but still compatible with 
the previous spectrum) and hints for the presence of a 
third peak around I ~ 850, which is not as high as ex- 
pected in a scenario with a large baryonic fraction. Si- 
multaneously the DASI experiment, which also found evi- 
dence for multiple peaks in the CMBR spectrum, gave an 
impressive and independent confirmation of a low baryon 
fraction, fl b h 2 = 0.022±g;gg| §, when sampling a differ- 
ent region of the sky and different frequencies. It is worth 
stressing that these high multipole data may still be af- 
fected by large systematic errors (see for example the 
consistency test in Table 3 in Ref. Q), thus all conclu- 
sions relying on them should still be taken with caution. 
This is especially true in view of the revised spectrum at 
£ > 300 from the Maxima-I experiment, which gives the 
wide range Q b h 2 = 0.0325±0.0125 @. 

Nevertheless it is important, on the basis of the new 
data now available, to undertake a detailed study of the 
compatibility of these data with SBBN. For this pur- 
pose we have performed, as in fl4| , a likelihood analysis 
of BOOMERanG/DASI CMBR data and SBBN in the 
parameter space (ft b h 2 , N u ), with N v the effective neu- 
trino degrees of freedom, and indeed we find a very good 
agreement. In particular the SBBN 95% CL region, cor- 
responding to N v = 2.8±0.3 and Q b h 2 = 0.020±0.004, 
has a large overlap with the analogous CMBR contour. 
This fact, if it will be confirmed by future experiments 
on CMBR anisotropy, can be seen as one of the greatest 
success, up to now, of the standard hot big bang model. 

As a byproduct of our analysis we also comment on 
the possible primordial 7 Li depletion, which has already 
been discussed in the literature BIIb&M . We find that 
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a depletion factor fa ~ 1/2 -j- 1/3 may reconcile observa- 
tions from Spite plateau with the value of flbh 2 . 

SBBN is well known to provide strong bounds on N v . 
On the other hand, Degenerate BBN (DBBN), first an- 
alyzed in Ref. p5|~p8| , gives very weak constraint on the 
effective number of massless neutrinos, since an increase 
in N v can be compensated by a change in both the chem- 
ical potential of the electron neutrino, fj, Ue = £, e T u , and 
flbh 2 . However, combining this scenario with the bound 
on baryonic and radiation densities allowed by CMBR 
data, it is possible to obtain rather strong constraints 
on N v even for DBBN. From our analysis we get the 
bound N v < 7, at 95% CL, when including Supernovae 
la (SNIa) data, which translates into a new and more 
stringent bound on background neutrino chemical poten- 
tials. 

Some caution is naturally necessary when comparing 
the effective number of neutrino degrees of freedom from 
BBN and CMBR, since they may be related to different 
physics. In fact the energy density in relativistic species 
may change from the time of BBN (T ~ MeV) to last 
scattering (T ~ eV). Specifically, if a neutrino has a 
mass in the range eV < m < MeV, and decays into ster- 
ile particles, like other neutrinos, majorons etc., with life- 
time f(BBN) < t < t(CMBR), then the effective number 
of neutrinos at CMBR would be sensibly different than 
at BBN [^9|. However, this possibility does not look too 
natural any longer, in view of the recent experimental re- 
sults on neutrino oscillation |3(],|3l| , showing that all ac- 
tive neutrinos are likely to have masses smaller than eV. 
One could instead consider sterile neutrinos mixed with 
active ones, which could be produced in the early universe 
by scatterings and subsequently decay. However, for mix- 
ing angle large enough to thermalize sterile neutrinos 
p2| , one needs a sterile to active neutrino number den- 
sity ratio n s /nu ~ 4-10 4 sin 2 29 (m/keV)(10.75/. 9 *) 3/2 of 
order unity [|33| (9 is the mixing angle, and g* is the 
number of relativistic degrees of freedom). Hence us- 
ing the decay time, r w 10 20 (keV/m) 5 / sin 2 29 sec, one 
finds t ~ 10 17 (keV/m) 4 yr, which is much longer than 
the age of the Universe, so they would certainly not have 
decayed at t(CMBR). Seemingly a sterile neutrino with 
mass of few MeV would have the right decay time, but 
this is excluded by standard BBN considerations p4| , |35" 
Let us emphasize that even though the simplest mode 
allow to directly combine BBN and CMBR results, nev- 
ertheless one may consider more exotic scenarios |5(||l8| , 
where £lbh 2 changes between BBN and CMBR epochs, 
or quintessence, which would result in a change of N v 
between BBN and CMBR |f). 

The paper is organized as follows: Section II is devoted 
to a brief review of the data used in our analysis, which 
is contained in Section III. Finally in Section IV we give 
our conclusions. 



II. BBN AND CMBR DATA 

A faithful estimate of primordial Deuterium is pro- 
vided by Ly-a features in several Quasar Absorption 
Systems (QAS) at high red-shift [z > 2). The most 
recent analysis of a four QAS sample gives D/H = 
(3.0±0.4)xl0~ 5 [3S[] . A new measurement has been also 
presented from observations of the Q-2206-199 QAS, at 
red-shift z ~ 2, which gives D/H = (1.65±0.35)xl(r 5 
p9| . When combined with the data of Ref. this re- 
sult gives a sensibly lower estimate for D abundance, 
D/H = (2.2±0.2)xlCT 5 . Nevertheless^ as reference 
value we will use the result quoted in ]3q] , but we will 
comment in our final discussion on the possible impact 
of what was found in |39| in the determination of Of,/i 2 
from BBN and CMBR new data. 

For the A He mass fraction, Yp, the key results come 
from the study of HII regions in Blue Compact Galaxies. 
The most complete and homogeneous sample has been 
analyzed in Ref. @, giving the value Y P = 0.244±0.002. 
A recent study, however, has pointed out the presence 
of possible systematic errors in inferring the total He 
abundance due to both imperfect ionization and non uni- 
form temperature distribution ]4lf| , leading to a typical 
overestimation of (2 — 4)% of Yp. This issue of course 
deserves a deeper study to understand if uncertainties in 
A He measurements are actually dominated by systematic 
effects. Notice that in the extreme value as low 

as Yp = 0.234 may represent a new problem for the very 
consistency of BBN scenario, in view of the low D re- 
sult of (3^1 . In what follows we will use with caution the 
result of Ref. ^(J quoted above. 

The estimate of 7 Li primordial abundance using Spite 
plateau can be spoiled by four possible systematic ef- 
fects 1 23 : a) Galactic Chemical Evolution (GCE), which 
is poorly known; b) corrections for possible depletion 
of initial star surface abundance; c) the very method 
of how 7 Li is obtained from Spite plateau; d) pres- 
ence of anomalous stars in the samples. In particular 
the effect due to GCE was long assumed to be negli- 
gible for metal poor stars in view of its apparent uni- 
formity, but this has recently been questioned due to 
observation of some amount of Be. Furthermore, data 
shows a statistically significant increase with Fe/H, as 
shown in p3| , leading to a primordial Lithium abundance 
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Evidence for this effect was 



instead missing in a previous analysis j42j, where it was 
found 7 Li/H = (1.73±0.21)xlCT 10 . The effects b) and 
c) have also recently been studied in Q, where it is 
pointed out that Spite plateau can be well reproduced 
by models with a strong diffusion effect, and would be a 
factor two lower than the primordial abundance. 

For these reasons, at present it is not appropriate to 
include 7 Li in a likelihood analysis of BBN. As in [|TJ, we 
will rather estimate from BBN prediction the depletion 
factor fa= 7 Liobs/ 7 Li pr i m , using as a reference result the 
one quoted in [E3|. 
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The anisotropy power spectrum from BOOMERanG 
experiment was estimated in 19 bins between I = 75 and 
I = 1025. Since the correlation matrix still is not public 
available, we will assume the data points to be indepen- 
dent. The data provide evidence for the presence of 3 
peaks at £ ~ 210±|, 550±^ 2 , 840±^ 3 , with an amplitude 
of ~ 72fiK, 49[iK and AbfxK respectively f43| . In our 
analysis we include a 10% correlation between the signal 
Cb in the bin B, a calibration uncertainty of 25% in AT 2 
and a gaussian uncertainty of 1.4' in the beam. Further- 
more, since the signal at very high multipoles [l > 850) 
could be severely affected by the presence of systematic 
effects, we apply a jackknife test repeating the analysis 
without these datapoints, finding no significant changes 
in our results. For the DASI data we include the win- 
dow functions available on the corresponding web site 
pij . We also include a 8% calibration error. There is an 
~ 20% overlap of the two regions of the sky covered by 
the two experiments but we do not take this effect into 
account in our analysis. In fact we believe that this corre- 
lation should not affect our conclusions, since our result 
appears stable when removing the DASI data points. 

III. LIKELIHOOD ANALYSIS 

The likelihood analysis of the BBN data has been per- 
formed using the method already described in details 
in |f(|. To constrain the values of the parameter set 
(N u ,n b h 2 ), for SBBN, and (£ e , N v , il b h 2 ), for the degen- 
erate scenario, from the data on 4 He and D, we define 
the likelihood function Cbbn = Ljj L<L He , where each 
likelihood function, assuming gaussian distribution for 
the errors, is given by the overlap of a theoretical and 
experimental distribution, 

The Yf x and af x are the experimental results and 1 — a 
errors for the i-th nuclide, Yf 1 the theoretical predictions 
obtained by an updated BBN code developed over the 
last few years p5| , p!6[ . Finally, the theoretical af 1 can 
be found by linear propagation of the uncertainties of 
the various nuclear rates entering in the nucleosynthesis 
reaction network |p2fl . 

For the BOOMERanG and DASI experiments we ap- 
proximated the likelihood function of the CMBR signal 
inside the bins, Cb, as a gaussian variable. The likeli- 
hood for a given cosmological model is then defined by 
-2\hC C mbr = (Cf - C B x f/a%, where Cf is the the- 
oretical signal. Our database of models is sampled as in 

s can be seen from Figure 1, the dotted (red) line, 
which represents the 95% CL contour of SBBN, is in very 



good agreement with new CMBR data, and the Sl b h 2 
tension between primordial nucleosynthesis and CMBR 
anisotropy seems to be completely solved. The constraint 
on fl b h 2 (N u ) can be obtained by marginalizing the total 
likelihood function C — Csbbn'£>cmbr with respect to 
N u (il b h 2 ). By this procedure we get the two estimates 
fl b h 2 = 0.019±0.003 and N v = 2.8±0.4, both at 95%. 

The result on N v beautifully suggests the simplest 
scenario of three light active neutrinos. It is there- 
fore perfectly meaningful to fix from the very beginning 
N v = 3.034 (Mj-gjj], which leads to the same interval for 
il b h 2 . In particular, for Q b h 2 = 0.019 the nuclei abun- 
dances evaluate to D/H = 3.26xl0" 5 , Y P = 0.2471 and 
7 Li/H = 3.31xl0~ 10 . 




0.010 0.015 0.020 0.025 0.030 0.035 0.040 

n b h 2 

FIG. 1. The 95% CL contours for degenerate BBN 
(dot-dashed (red) line), new CMBR results only with age 
prior, t > llgyr (full (green) line), and only with SNIa prior 
(dashed (blue) line) are shown. The combined analyses cor- 
respond to filled areas: DBBN + CMBR + age (light (green) 
region), DBBN + CMBR + SNIa (dark (blue) region). The 
dotted (red) line is the 95% CL contour of SBBN. 



Although SBBN seems to be in very good agreement 
with the new CMBR data, providing a very strong con- 
straint on the neutrino degrees of freedom, it relies on 
the theoretical assumption that background neutrinos 
have negligible chemical potential, as their charged lep- 
ton partners. Even though this hypothesis is perfectly 
justified by Occam razor, models have been proposed in 
the literature p9| , |50| -p6[ where large neutrino chemical 
potentials can be generated. It is therefore an interest- 
ing issue for cosmology, as well as for our understand- 
ing of fundamental interactions, to try to constraint the 
neutrino-antincutrino asymmetry from cosmological ob- 
servables. It is well known that degenerate BBN gives 
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severe constraints on the electron neutrino chemical po- 
tential, —0.06 < £ e < 1.1, and weaker bounds on the ones 
of both a and r neutrino, |£(U,t| < 5.6 -j- 6.9 p^j. This 
occurs since electron neutrinos are directly involved in 
neutron to proton conversion processes which eventually 
fix the total amount of 4 He produced in nucleosynthesis, 
while only enters via their contribution to the expan- 
sion rate of the universe. Combining this scenario with 
the bound on baryonic and radiation densities allowed 
by CMBR data, it is possible to obtain rather stronger 
constraints on all these parameters. Such an analysis 
was previously performed in |14|,|57| using BOOMERanG 
and MAXIMA data of Refs. We recall that neu- 

trino chemical potentials contribute to the total neutrino 
effective degrees of freedom N v as 
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with 8 V the contribution of relativistic degrees of freedom 
other than neutrinos and photons. Notice that, in order 
to get the most stringent bound on £ a we have to as- 
sume that all relativistic degrees of freedom, other than 
photons, are given by three active (possibly) degenerate 
massless neutrinos, i.e. 8 V = 0. Similarly the upper limit 
on 8 V can be obtained from the results of our analysis in 
the case ^ M , r = 0. We stress that in any case a value for 
N v sensibly different than three does require a non van- 
ishing chemical potential for electron neutrinos, or more 
generally a non thermal spectrum. 

Figure 1 summarizes our main results for the DBBN 
scenario. Defining AN V — N v — 3, we plot in the plane 
(AN V , fl b h 2 ) the 95% CL contour allowed by DBBN (dot- 
dashed (red) line), together with the analogous 95% CL 
region coming from the CMBR data analysis, with only 
weak age prior, r > llgyr (full (green) line). Finally, 
the light (green) filled region is the 95% CL region of the 
joint product distribution C = Cdbbn'£-cmbr- The 
main new feature, with respect to the results of Ref. JlJ] 
is that the resolution of the third peak shifts the CMBR 
likelihood contour towards smaller values for r2{,/i 2 , so, 
when combined with DBBN results, it singles out smaller 
values for N v . In fact from our analysis we get the bound 
N v < 8, at 95% CL, which translates into the new bounds 
-0.01 < £e < 0.25, and, for 8 V = 0, |^ iT | < 2.9, sensi- 
bly more stringent than what can be found from DBBN 
alone. 

A similar analysis can be also performed combining 
CMBR and DBBN data with the Supernova la data Q , 
which strongly reduce the degeneracy between f2 m and 
fJ A - At 95% CL we find (dark (blue) filled region in Fig- 
ure 1) N v < 7, corresponding to —0.01 < £ e < 0.22 and 
|£/i,r| < 2.6, for 8 V = 0. In the other extreme scenario, 
where basically all extra contributions to Hubble param- 
eter are given by extra relativistic species, we get 8 V < 4. 

Another possibility to break the degeneracy between 
O m and Cl\, is to put priors on the age of the uni- 
verse t, as pointed out in Ref. [^tJ. In Figure 2 we 



show the normalized likelihood functions for age priors of 
t > 10, 11, .., 14 gyr, using only CMBR data. It is clear 
that one needs the slightly unrealistic prior of r > 13 gyr 
to get bounds stronger than AN < 4, as obtained by the 
inclusion of SNIa data. 



T 
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FIG. 2. CMBR likelihoods for age priors of r > 10, 11, .., 14 

gyr- 



Recently Ref. |39| stressed the point that the inclusion 
of large scale structure data can provide a lower bound on 
AN V . In Ref. |)^] SNIa data are not considered, neither 
a DBBN scenario; however, where comparison is possi- 
ble, the results on upper bounds are in fair agreement 
with our results]]. It is worth noticing once again that 
our rather stringent bound on N u is the outcome of the 
combined analysis of CMBR and DBBN. Each of the 
two corresponding likelihood contours in fact, taken sep- 
arately, give a much weaker bound (see Figure 1). 

As we mentioned in Section 2, it has recently been 
stressed that depletion effects on 7 Li may be efficient 
in reducing the primordial abundance down to the value 
observed in Spite plateau. In Figure 3 we plot the 7 Li de- 
pletion factor pjr, defined as the ratio of the experimental 
value of Ref. |2!| and the theoretical estimate from our 
BBN code. Values for fa of the order of 1/2 -j- 1/3 cannot 
be ascribed to a statistical fluctuation in the star sample 



*A study of the effects of large neutrino asymmetries on 
CMBR + large scale structures has also been performed in 
Ref. @, where the compatibility with DBBN for ra„ < leV 
has been analyzed. Their results refer to a critical universe 
with no cosmological constant. 



4 



considered in |£3| , but should rather be understood by a 
careful analysis of all systematic effects which we briefly 
reviewed in Section II. 

There are some points we would like to address as final 
remarks. First of all we stress once again that further 
data on the third peak in the CMBR anisotropy spectrum 
are needed to check for possible systematics. This is a 
crucial point for a clean determination of the baryonic 
fraction, since discrimination between SBBN and DBBN, 
or SBBN and other theoretical framework for light nuclei 
production, relies on both second and third peak heights. 
In this respect we note the good agreement between the 
BOOMERanG and DASI results. 

As a second observation, we recall that we already 
pointed out that a new measurement of primordial D 
has been reported recently, leading to a weighted aver- 
age D/H = (2.2±0.2)xl(T 5 §§. If we adopt this dif- 
ferent estimate, the overlap of SBBN and CMBR con- 
tours decreases, but still there is a good agreement at 
95% CL. Using the SBBN likelihood analysis only we 
find in this case, at 95% CL, Vl b h 2 = 0.0241°;°°^ and 
N v = 2.7±0.3, while combining this result with' CMBR 
data and using the joint likelihood distribution C we get 
n b h 2 = 0.023±g;gg| and N v = 2.7±0.3. 




0.016 0.018 0.02 0.022 0.024 



ft B h 2 

FIG. 3. The 7 Li depletion factor, defined as the ratio be- 
tween the experimental and theoretical values. 



of information for the whole microscopic theory of funda- 
mental interactions. The increasing precision in measure- 
ments of primordial abundances of light nuclei, and the 
impressive progress in measuring the CMBR anisotropy, 
are conspiring to give us a very precise determination 
of N v . Despite of the conservative expectation of three, 
light, active neutrinos, largely non-degenerate, it should 
be stressed that many other scenarios have been consid- 
ered in the literature, based on theoretical ideas which, 
going beyond the Standard Model, try to grasp possible 
extension of our knowledge of fundamental interactions 
at higher energy. It is really exciting that, along with 
customary accelerator physics, we have at hand a severe 
way to scrutinize these models by cosmological measure- 
ments. 

In this paper we have studied in details the implica- 
tions of the new BOOMERanG and DASI data on CMBR 
anisotropy for the estimation of the baryonic energy den- 
sity fraction, compared with the predictions of standard 
BBN, in the parameter space Vl b h 2 — N v . Observation 
of the third peak at multipole I ~ 850 turned into a 
sensible improvement of the compatibility of the two in- 
dependent ways of constraining fl b , and single out the 
values fl b h 2 = 0.019±0.003, and N v = 2.8±0.4, both at 
2cr. 

We have also considered the scenario of a degenerate 
neutrino background, which strongly affects primordial 
nuclei production. The new CMBR BOOMERanG and 
DASI data lead to a new and stronger constraint on the 
effective relativistic degrees of freedom, N v < 8 (only 
weak age prior), or N v < 7 (with only SNIa prior), both 
at 95% CL, which bounds more severely the neutrino 
chemical potentials, |£ M , T | < 2.9, and — 0.01 < £ e < 0.25, 
and I^tI < 2.6, and —0.01 < £ e < 0.22, respectively. 
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IV. CONCLUSIONS 

It is a great success of cosmology and astrophysical ob- 
servations that severe constraints can be put on the num- 
ber of neutrino degrees of freedom and, more generally, of 
light particle species which were relativistic at the epoch 
of recombination. Of course this is a fundamental piece 
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